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COMPUTER AIDED SOLUTION OF THE INVARIANCE 
EQUATION FOR TWO- VARIABLE GINI MEANS 

SZABOLCS BAJAK AND ZSOLT PALES 

Abstract. Our aim is to solve the so-called invariance equation in the class of two- 
variable Gini means {G P ,g : p,q £ R}, i.e., to find necessary and sufficient conditions 
on the 6 parameters a, b, c, d, p, q such that the identity 

Gp,q(G a ,b{x,y),G Ct d(x,y)) = G Piq (x,y) (x,y G K+) 

be valid. We recall that, for p ^ q, the Gini mean G P:q is defined by 

The proof uses the computer algebra system Maple V Release 9 to compute a Taylor 
expansion up to 12th order, which enables us to describe all the cases of the equality. 



1. Introduction 

Let M + denote the set of positive real numbers throughout this paper. A two- variable 
continuous function M : — > M + is called a mean on M + if 

(1) min(x, y) < M(x, y) < max(x, y) (x, y G M + ) 

holds. If both inequalities in ([T|) are strict whenever x ^ y, then M is called a strict 
mean on M + . 

Given two means M, iV : — > M. + and x, y G M + , the iteration sequence 

x x :=x, V\'=V, 
x n+1 := M(x n , y n ), y n+1 := N(x n , y n ) (n G N) 

is said to be the Gauss-iteration determined by the pair (M, N) with the initial values 
(x,y) G Ml. It is well-known (cf. [T], [2]) that if M and are strict means then the 
sequences (x n ) and (y n ) are convergent and have equal limits M ® A r (x, which is a 
strict mean of the values x and y. The mean M®N defined by this procedure is called 
the Gauss composition of M and N. 

A key result characterizing the Gauss composition of means is the following state- 
ment: Given two strict means M, N : IR^ — > R + , their Gauss composition K = M ® N 
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is the unique strict mean solution K of the functional equation 

(2) K(M{x,y),N{x,y))=K{x,y) (x,y Gl + ) 

which is called the invariance equation. 

The simplest example when the invariance equation holds is the well-known identity 



x + y 2xy 

x y = \^> r~ M eR +)> 

V 2 x + y 
that is, 

5(x, y) = 9 (./l(x, y), J£(x, y)) (x, y G R + ), 

where .A, 9, and !K stand for the two- variable arithmetic, geometric, and harmonic 
means, respectively. Another less trivial invariance equation is the identity 

A ® S(x,y) = A® $(A(x,y), $(x,y)) (x,y G R+), 

where A <g> 9 denotes Gauss' arithmetic-geometric mean defined by 

A®9(x,y) = (- [ 2 dt ===] (x,yeR + ). 

x 71 Jo y x cos t + y sin £/ 

The invariance equation in more general classes of means has recently been studied 
extensively by many authors in various papers. The invariance of the arithmetic mean 
A with respect to two quasi-arithmetic means was first investigated by Matkowski [3] 
under twice continuous differentiability assumptions concerning the generating func- 
tions of the quasi-arithmetic means. These regularity assumptions were weakened step- 
by-step by Daroczy, Maksa, and Pales in the papers |3], [5], and finally this problem 
was completely solved assuming only continuity of the unknown functions involved [2J. 
The invariance equation involving three weighted quasi-arithmetic means was studied 
by Burai |6j, [7J and Jarczyk-Matkowski [8J, Jarczyk |9j. The final answer (where no 
additional regularity assumptions are required) has been obtained in [9]. In a recent 
paper, we have studied the invariance of the arithmetic mean with respect to two so- 
called generalized quasi- arithmetic mean under four times continuous differentiability 
assumptions [10J. The invariance of the arithmetic mean with respect to Lagrangian 
means was the subject of investigation of the paper [llj by Matkowski. The invari- 
ance of the arithmetic, geometric, and harmonic means with respect to the so-called 
Beckenbach-Gini means was studied by Matkowski in [TJ] . Pairs of Stolarsky means for 
which the geometric mean is invariant were determined by Blasinska-Lesk-Glazowska- 
Matkowski [13J. The invariance of the arithmetic mean with respect to further means 
was studied by Glazowska- Jarczyk-Matkowski [T3] and Domsta-Matkowski |15| . 

An important class of two variable homogeneous means are the so-called Gini means 
(cf. Gini [16J). Given two parameters p, q G R, the two-variable mean G p ^ q : R^ — > R + 
is defined by the following formula 



G p Jx,y) 



X P + yP 

xi + yi 



i 

p-q 



for p ^ q, 



f x p \nx + y p \ny\ 

ex P V^~v for p = q, 

\ x p + y p J 

for x, y G R+. 

The class of Gini means is a generalization of the class of power means, since taking 
q = 0, we immediately get the power (or Holder) mean of exponent p. 
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The aim of this paper is to solve the invariance equation in the class of Gini means, 
i.e., to solve (J2J) when each of the means M, N, and K is a Gini mean. More precisely, 
we want to describe the set of all 6-tuples (a, b, c, d,p, q) such that the identity 

(3) G Ptq (G aib (x, y), G c4 (x, y)) = G m (x, y) (x, y G R+) 
holds. The main result of this paper is contained in the following theorem. 

Theorem. Let a,b,c,d,p,q G R. Then the invariance equation ([3]) is satisfied if and 
only if one of the following possibilities hold: 

(i)a + b = c + d = p + q = 0, i.e., all the three means are equal to the geometric 
mean, 

(ii) {a, b} = {c,d} = {p, q}, i.e., all the three means are equal to each other, 
(Hi) {a, b} = { — c, — d} and p + q = 0, i.e., G p<q is the geometric mean and G a ^ = 
G-c-d, 

(iv) there exist u, v G R such that {a, b} = {u + v,v}, {c, d} = {u — v,—v}, and 
{P) <?} = { M >0} (i n this case, G P)Q is a power mean), 

(v) there exists w G R such that {a, b} = {3w,w}, c + d = 0, and {p, q} = {2w,0} 
(in this case, G p<q is a power mean and G c ^d is the geometric mean), 

(vi) there exists w G R such that a + b = 0, {c,d} = {3w,w}, and {p,q} = {2w,0} 
(in this case, G PA is a power mean and G a ^ is the geometric mean). 

As an obvious consequence, we obtain the following solution for the so-called Matkowski- 
Suto equation, i.e., when G PA is equal to the arithmetic mean in ([3]). 

Corollary. Let a, b,c,d G R. Then the Matkowski-Suto equation 
G atb (x,y) + G Cid (x,y) = x + y (x,yeR + ) 

is satisfied if and only if one of the following possibilities hold: 

(i) {a, b} = {c, d} = {1,0}, i.e., the two means are equal to the arithmetic mean, 
(ii) there exist uel such that {a, b} — {1 + v, v}, {c, d} = {1 — v, —v}, 
(Hi) {a, b} = {§,5} an d c + d = (in this case, G Ci d is the geometric mean), 
(iv) a + b = and {c, d} = {|, |} (in this case, G a ^ is the geometric mean). 

In view of the homogeneity of Gini means, identity ([3]) is equivalent to the equation 

(4) g p , 9 (ga,t(e a ,e- a! ),G c , d (e a ,e- 3! )) = ^ 

G Pt q (g ; G ) 

The main and simple idea of the proof of our Theorem is to compute the Taylor 
expansion of the function F : R — > R defined by 

G m (G atb (e x ,e- x ),G c>d (e x ,e- x )) 



(5) F(x) 



Gp,q{ 



at the point x — 0. Using the symmetry of Gini means, it also follows that F is an 
even function. Therefore all coefficients Ck of odd order in the Taylor expansion of F 
are equal to zero automatically. The validity of fll]) yields that all of the even order 
coefficients also vanish. To derive the necessity of the conditions in our Theorem, 
we will only use the equalities C2 = ■ ■ ■ = Cyi = 0, which produce 6 equations for 
the 6 unknown parameters a, b, c, d, p, q. The sufficiency of the conditions so obtained 
will be proved by a simple argument which, implicitely, shows that the equalities 
C2 = ■ ■ ■ = C12 = yield C^fc = for k > 7. Unfortunately, these 6 equations are 
so complicated that their solution requires the power of the Maple computer algebra 
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package. This will be done in the next section, where we also show the sequence of 
Maple commands that were performed during our computation. If the interested reader 
executes all these commands in a Maple worksheet, then all the computations can be 
repeated and checked again. During this computation, we analyze each of the Taylor 
coefficients separately and use the information so obtained at the subsequent steps. 



2. The proof of Theorem 
First we recall the characterization of the equality of two variable Gini means. 
Lemma. ( Cf. |17| ) Let a, b, c, d G R. Then the identity 

(6) G a>b (x, y) = G Cid (x, y) (x, y G R+) 

holds if and only if one of the folowing possibilities is valid: 

(i) a + b = c + d = and, in this case, the two means are equal to the geometric 
mean, 

(ii) {a, b} = {c, d}. 

Proof of the Theorem. Assume that ([3]) holds. Then the function F defined by (jSJ) is 
identically 1 on R. Thus, for the A;th-order Taylor coefficient Ck defined below, we 
have: 

(7) C k := = (Jfe G N). 

Since F is even, thus Ck = for all odd k G N. 

In view of the symmetry of the Gini means in the parameters, we may assume that 
a > b, c > d, and p > q in the sequel. 

For the calculations in Maple, define the Gini mean G p ^ q (x, y) = G(p, q, x, y) and the 
function F given by by performing the following commands: 

> G:=(p,q,x,y)->((x~p+y~p)/(x~q+y~q))~(l/(p-q)) ; 

F:=x->(G(p,q,G(a,b,exp(x) ,exp(-x)) ,G(c,d,exp(x) ,exp(-x)))) 
/ (G(p,q,exp(x) ,exp(-x))) ; 

This yields the following output: 

' x p + y % 



G : = (p, q, x, y) \-t 
F : = x \-> 



x q _|_ y q ^ 

G(p, q, G{a, 6, e x , e~ x ), G(c, d, e x , e~ x )) 



G(p,g,e x ,e x ) 

In the first step, we evaluate the 2nd-order Taylor coefficient C<i by executing 
> C[2] :=simplify(coeftayl(F(x) ,x=0,2)) ; 
This results 

1 k 1 1. 1 1 
C 2 :=-a + -b + -c+- A d--p--q 

Therefore, by C% = 0, we obtain our first necessary condition: (a + b + c + d)/4 = 
(p+q)/2. If one tries to compute C4, Cq, . . . , then the expressions obtained are so com- 
plicated that it is hard to get further information. In order to simplify the evaluation 
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of the higher-order Taylor coefficients, we introduce the following notations 

a + b + c + d p + q 



w 



4 

a + b — (c + d) 
' 4 ' 

2 

p-q 



r :- 



t :-- , 

2 

(a-bf + (c-d) 2 
8 

__ {a-b) 2 - (c-rf) 2 
S '~ 8 

(In the definition of w we utilized the condition Ci = 0.) Then we can express the 
parameters a, b, c, d,p, q in the following form: 

> a: =w+v+sqrt (r+s) ; b : =w+v-sqrt (r+s) ; 
c : =w-v+sqrt (r-s) ; d : =w-v-sqrt (r-s) ; 
p:=w+sqrt(t) ; q: =w-sqrt (t) ; 



W + V + y/r + s 



w + v — y/r + s 



w — v + y/r — s 



w — v — y/r — s 
p := w + yft 
q := w — yft 

Now we evaluate the 4th order Taylor coefficient by inputting: 

> C[4] :=simplify(coeftayl(F(x) ,x=0,4)) ; 
Then we obtain 

La := -tw vs wr 

3 3 3 

The condition C4 = yields that wt = wr + vs. 

If w = 0, then p + q = and hence G PA is equal to the geometric mean. Therefore, 
the invariance equation can be rewritten as 

G atb (x,y)G Cjd (x,y) = xy (x,y eR+). 

This results 

G at b(x, y) = — - = G^ d {x, y) (x, y G 

G c , d {l/x, l/y) 

Using the Lemma again, this identity yields that either a + b = c + d = 0oT {a, b} = 
{— c, —d} must hold. Together with p + q = 0, these equations show that either con- 
dition (i) or condition (iii) of our theorem must be satisfied. Conversely, if conditions 
(i) or (iii) hold then an easy computation yields that ([3]) is satisfied. 

In the rest of the proof, we assume that w is not zero. Then, we can express t in 
terms of w, v, r, s: 

> t:=r+v*s/w; 



6 SZ. BAJAK AND ZS. PALES 

, x vs 
8 t:=r + — 

w 

Next, we evaluate the 6th order Taylor coefficient: 

> C[6] :=simplify(coeftayl(F(x) ,x=0,6)) ; 

— 2(— 3w 2 s 2 + 3v 2 s 2 — 15w 2 rv 2 — 5w 3 sv — 10v 3 sw + 15w 4 v 2 ) 

Ok '.— 

A5w 

If v = 0, then the condition Cq = simplifies to w 2 s 2 = 0, whence s = follows. 
Therefore (JBJ) yields t = r and we obtain that a = c = p and b = d = q which means 
that condition (ii) of our theorem must be fulfilled. In this case, it is obvious that ([3]) 
is satisfied. 

In the rest of the proof, we assume that v is also not zero. Observe that the 6th 
order coefficient Cq does not involve higher-order powers of r. Therefore, the equation 
Cq = can be solved for r. Temporarily, we denote this solution by R: 

> R:= (15*w~4*v~2-3*w~2*s~2+3*v~2*s~2-5*w"3*v*s-10*w*v~3*s)/(15*w"2*v"2) ; 

/ n \ 15w A v 2 — 3w 2 s 2 + 3v 2 s 2 — 5w 3 vs — 10wv 3 s 

W R := u 2 2 

low 

Finally, we evaluate the 13th order Taylor polynomial of F(x) at x = (the Maple 
output is suppressed by putting : instead of ; to the end of the Maple command, for 
the sake of brevity), then we extract the 8th, 10th and 12th order Taylor coefficients, 
replace r by R and factorize the expressions so obtained by inputting: 

> T:=simplify(taylor(F(x) ,x=0,13)) : 
for i from 8 to 12 by 2 

do C [i] : =simplif y(subs (r=R, simplif y(coef f (T,x, i) )) ,f actor) od; 

„ (V — W)(V + W)S / o A r q o a O o 

C 8 := K — I K — ' (2100wV - 3850u>Vs + 4200u>V - 255wv 3 s 2 + 153v 2 s 3 
70875urt> 2 

- 9245wVs - 7395w V - 153wV) 

„ 2(v — W)(V + W)S / r Q A o Q 79 77 

C 10 := -i ^ — — ^(28500wV - 59675wVs + 34470wVs 2 + 20100wV 

1063125w 5 f 4 

- 299575w% 6 s - 4260wVs 3 - 73200wVs 2 - 930«wV + 66600wV 
+ 4805wVs 3 - 286500wVs + 279vV - 169020wVs 2 - 16740wVs 4 

- 558w Vs 5 + 45955wVs 3 + 17670w V + 279wV) 

C 12 := 2{V ~ W) }l + ^ (- 3272692500^%^ + 2218110(WW - 54365475?/; W 
2631234375w'f b 

- 25317375wVs 3 + 335826«A>V - 559710Ws 6 - 215221875wV 2 s 

+ 22875570wVs 5 + 16977870wVs 6 - 7649370wVs 6 - 1246797750wV°s 

- 34684335wVs 5 - 777170000W 1 Vs 2 - 159926550wVs 4 - 335826wVs 7 

+ 641072375w 1( W + 270963000wW - 1046615000wW + 11659365wVs 5 

- 133190250wVs 4 - 1967022000w 1 Vs + 177650700w Vs 4 - 8768790wW 
+ 385915750wW + 149400wVs 5 - 98002025m; W + 76725000w 7 t> 13 

- 57172500U7 9 /; 11 - 478665000w 1 V + 188100000m; 13 i; 7 - 111942w 6 s 7 + 111942«V) 



INVARIANCE EQUATION FOR TWO VARIABLE GINI MEANS 7 

(In the entire Maple computation, this is the only step which requires a considerable 
processing time. On a DualCore 2.6 GHZ processor, the elapsed time was less than 1 
minute.) 

Cg, Cio and Cn are obviously zero if s(v — w)(v + w) = 0. Thus, we have to consider 
the three subcases: s — 0, v — w, and v = —w. 

In the case s = 0, flSj) and Q imply that t = r = w 2 . Then we get that 

{a, b} = {2w + v,v}, {c, d} = {2w - v, —v}, {p, q] = {2u>, 0}, 

i.e., condition (iv) holds with u := 2w. Conversely, if condition (iv) holds and u ^ 0, 
then we have 

G P}q (G ajb (x,y),G Cid (x,y)) = G ufi (G u+v<v (x,y),G u - v - v (x,y)) 

( x u+v _|_ yU+v x u ~ v + y u ~ v \ » / x u+v + y uJrV x u y v + y u x v \ " 
2(x v + y v ) + 2(x- v + y- v )J ~ \ 2(x v + y v ) + 2{x v + y v ) ) 

— = — J = G ufi (x,y) = G PyQ (x,y). 

Thus ([3]) is satisfied if u ^ 0. If u = 0, then the parameters also fulfil condition (iii), 
hence ([3]) holds in this case, too. 

If v = w, then, by (jSJ) and (Q, r = w 2 — s and t = r + s = w 2 , respectively. Hence, 

{a, b} = {3w, w}, c + d = 0, {p, q} = {2w, 0}, 

i.e., condition (v) holds. Conversely, if condition (v) holds and w^0, then, using the 
identity (1101) with u := 2w, v := w, we have 

G Pjq [G a! b(x,y),G Ct d(x,y)) = G2 W! o(G2 W +w,w(x,y),G2 W - w - w (x,y)) 

= G 2w ,o{x,y) = G Piq (x,y), 

which shows that is fulfilled. If w — 0, then all the three means are geometric 
means, and hence ([3]) holds trivially. 

The last case when v = —w holds, similarly to the case v = w, implies that condition 
(vi) is valid. If condition (vi) holds, then ([3]) can also be verified. 

In the rest of the proof, we can assume that s(v + w)(v — w) is not zero. The Taylor 
coefficients Cg, Cio, and G\i are of the form 

_ (v — w)(v + w)s _ 2(v - w)(v + w)s _ 2(v - w)(v + w)s 

8 ~ 70875^ 8 ' 10 _ 1063125w 5 i; 4 10 ' 12 ~ 2631234375w 7 w 6 15 

where P 8 , P 10 , and Pyi are polynomials of the variables v,w,s. They can be obtained 
by the following Maple commands (whose output is suppressed): 

> P[8] :=op(5,C[8]) : P [10] : =op(5 , C [10] ) : P [12] :=op(5,C [12] ) : 

The equalities C$ = C w = C\i = and s(v + w)(v — w) ^ imply that P s = P w = 
Pi 2 = 0. In what follows, we show that there is no solution v,w,s to this system of 
equations. 

The variable s is a common root of the polynomials Pg and Pio- Therefore the 
resultant Pg,io of these two polynomials (with respect to s) is zero: 

> R[8,10] :=f actor (resultant (P [8] ,P[10] ,s)); 
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# 8i io :=136687500?/V 5 (t; - w) 2 (v + wf 

(1178440166794705680 v 18 - 34849488132334981400 w 2 v 16 

+ 27095657773476976150 u>V 4 + 2157163953185024831539 v 12 w 6 

+ 19335728720363587723895 w 8 v 10 + 77098340762854904758838 w 10 v 8 

+ 135541716064734053550290 w 12 v 6 + 52974528518488497499557 v 4 w u 

+ 2100034048587009260985 w w v 2 + 44498612407766474466 w 18 ) 

Since vw(v — w)(v + w) ^ holds, we get that v and w are solutions of a homogeneous 
two variable polynomial equation of degree 18. Writing w in the form 

> w:=z*v; 

w := zv 

we get that z is a root of a 18th degree polynomial Ps,io (more usefully, z 2 is a root of 
a 9th degree polynomial), where: 

> P[8,10] :=simplify(op(4,R[8,10])/v~18); 

P 8 ,io :=1 178440166794705680- 34849488132334981400/ 

+ 27095657773476976150/ + 2157163953185024831539/ 
+ 19335728720363587723895/ + 77098340762854904758838/ 
+ 135541716064734053550290/ 2 + 52974528518488497499557/ 4 
+ 2100034048587009260985/ 6 + 44498612407766474466^ 18 . 

The variable s is also a common root of the polynomials P 8 and P 12 . Therefore the 
resultant i?8,i2 of these two polynomials (with respect to s) is again zero: 

> R [8 , 12] : =f actor (resultant (P [8] , P [12] , s) ) : 

Since vw(v — w)(v + w) ^ holds, we now get that v and w are solutions of a 
homogeneous two variable polynomial of degree 26, whence we get that z is a root of 
the 26th degree polynomial Ps,i2) where: 

> P[8,12] :=simplify(op(4,R[8,12] )/v~26) ; 

P[8, 12] := 8196063700595383871701091232 - 179090512353635410423157248720/ 
- 2262574745604112043731392907114/ + 11198535065282946302316347517923/ 
+ 369075355861065090753396085824722/ + 3321203212966063219800014204539694/ 
+ 17018221168597358591328346358640128/ 2 + 55161742271395394206883716537690208/ 4 
+ 113024609788553283598449985201081964/ 6 + 136472191224999845881431378284988722/ 8 
+ 83840233563357841801204648333566258/ + 19391722782753178903737004919064981/ 2 
+ 1234978033803167388960240130106010z 24 + 95711050739605210548400442203992/ 6 

Now computing the resultant of the two polynomials Ps,io and Ps,i2 by 

> Q : =resultant (P [8 , 10] , P [8 , 12] , z) : 

it follows that Q is a (huge) nonzero number, hence Ps^o and Ps,i2 cannot have a 
common root. This proves that Ps.io and Ps,i2 can be simultaneously zero if and only 
if vw(v — w)(v + w) =0 holds. Hence C 8 = C w = C± 2 = can only hold also in this 
case. □ 
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